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Abstract 


This paper focuses on the result of inclined angle on bioconvection 
of porous media bounded by cavity wall square enclosure filled with 
both nanofluid and gyrotactic microorganisms passing through the media 
with pores. The dimensionless velocity, temperature, concentration, and 
mass transformation equations are solved by using the weighted residual 
Galerkin’s finite element method. The result of the inclination angle from 
56 = 0° to 6 = 180° in a square cavity is interpreted. The outcomes of incli- 
nation on various key parameters, such as Rayleigh number, bioconvective 
Rayleigh number, Peclet number, Brownian motion, and the ratio of buoy- 
ancy, are discussed. Furthermore, the mean Nusselt number, Sherwood 
number, and density number are discussed at vertical walls. 
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1 Introduction 


The present problem is dealing with the nanofluid flow through the porous 
square cavity with temperature difference, which has a wide range of applica- 
tions in recent years, such as geophysics, geothermal energy utilization, and 
many technologies. The bioconvection of the nanofluid containing gyrotactic 
microorganisms has a wide range of practical applications, such as chemical 
catalytic converters, buried electronic cables, pollutant dispersion in aquifers, 
food industrial forms, and so on. These types of many areas of applications 
are documented in these references [10, 24, 14, 15, 31]. 

The properties and utilization of the nanofluid were first introduced by 
Choi and Eastman [11] at ASME annual meeting. Many people have de- 
scribed the properties of nanofluids, such as [13, 30, 24, 23, 26]. In many 
electronic devices, like computers, boilers, converters, and so on, the angle of 
inclination to the surface affects the gravity force on the fluid, temperature 
gradient, and velocity of the fluid flow. In [32], the author expressed the free 
convection of the composite wall enclosure. Kuyper et al. [20] studied the 
effect of inclined angle on different flows in square cavity walls. Kuznetsov 
[21] explained the microscopic convection motion of the oxytactic microor- 
ganisms due to the temperature effect. Shermet and Pop [28] expressed that 
the result of the thermal movement of microorganisms in the nanofluid hav- 
ing gyrotactic microorganisms is a closed porous square cavity. Aziz, Khan, 
and Pop [5] presented the flow behavior of the nanofluid with gyrotactic mi- 
croorganisms on a flat plate. At viscous dissipation, the behavior of oxytactic 
microorganisms in porous square cavities was discussed [2, 3]. Jamuna and 
Balla [17] discussed the behavior of the heat source and sink of the gyrotac- 
tic microorganisms in the square cavity. The activation energy effect on the 
gyrotactic microorganisms was discussed in [18]. The influence of Soret and 
Dufour on free convection of the fluid flow in the inclined four-side closed 
walls was explained in [8]. The MHD double-diffusion in the porous square 
enclosure with radiation and chemical reaction and the outcome inclination 
of the porous square cavity filled with gyrotactic microorganisms with the 
heat transformation was discussed in [6, 7]. Nanofluid movement in an in- 
clined square cavity with gyrotactic microorganisms at MHD free convection 
was reported in [27]. 

The effect of angle movement of the square adiabatic wall on mixed con- 
vection of the nanofluid was explored in [16]. Aounallah et al. [4] explained 
the turbulent flow behavior of the nanofluid in an inclined square cavity on 
free convection, and Sheremet, Grosan, and Pop [27] investigated the free 
convective flow of nanofluid in inclined four-sided chamber with gyrotactic 
microorganisms. Tsai, Li, and Lin [29] discussed the inclination of the plate 
shield, and Aboueian-Jahromi, Hossein Nezhad, and Behzadmehr [1] studied 
the steady flow in inclined cylinders. 

Rajarathinam and Nithyadevi [25] examined the movement of Cu-water 
nanofluid in inclined cavity walls with pores. The thermosolutal Maragoni 
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effects of the bioconvective fluid flow with gyrotactic microorganisms on in- 
clined sheets were explained in [19]. Recently, Varol, Oztop, and Koca [34] 
explained different fluids’ laminar flow in the different inclined enclosures. 

Since, from the above literature survey, we note that many authors con- 
centrated on the inclined angle of different geometries with convection of 
nanofluid with gyrotactic microorganisms. The novelty of this paper con- 
tains the square-shaped cavity enclosure with fluid containing nanoparticles 
and gyrotactic microorganisms. Galerkin’s finite-element method is used to 
solve the nondimensional governing equations. 


2 Mathematical modeling 


We consider the bioconvection flow in an inclined two-dimensional porous 
four-sided square cavity of dimension L containing nanofluid with gyrotactic 
microorganisms. Let us assume that 6 is the inclination angle of the cavity 
wall with the horizontal surface. The vertical walls are maintained in various 
temperatures Tc and Ty, respectively (Ty>Tc ). The remaining walls were 
kept perfectly insulated. The direction of gravity force g acts opposite to the 
vertical axis (Y-axis). 

The steady-state Darcy—Boussinesq approximation governing equations 
are 


Ou Ov 


ey Ae ee Go 1 
oat Dy 0, (1) 
[Lb Op : 
pu = a, leper) (C—Cmin) —(1—-Cmin) p¢8(P—Te) + ynAp|gsin 4, (2) 
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Here y is the mean volume for microorganisms, Ap = pce — pr is the 
density difference of cell, T¢ is the cold wall temperature, Ty is the hot wall 
temperature, @,, is porous medium thermal diffusivity, C is the concentra- 
tion of nanoparticles, Co is nanoparticles average density, Cinin is minimum 
concentration of oxygen essential for microorganisms, C, is specific heat at 
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constant pressure, D,, is microorganisms diffusion coefficient, Dg is the Brow- 
nian diffusion constant, Dr is the thermophoretic diffusion coefficient, n is 
motile density number of microorganisms, g is the gravity force, chemotaxis 
constant is b, and the maximum speed of cell swims is wo. The average 
swimming velocities of microorganisms are u and v given as 


ae bwo OC Zs bwo 0c 
u= v= —- —. (7) 
AC dx’ AC Oy 
Consider dimensional stream function ~. Then u and v in x and y directions 


are considered as u = se and v = oe by introducing the boundedless vari- 


ables X= #,Y=4,0=~.6 Tar, Te, 6 = S Sno, and N = 
ng is the microorganism averaged denciey 
Substituting above unbounded variables into equation (1)—(7), then we 


get the following partial differential equations: 


where 
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We have WV = 0 for all sides, 
g¢=1,0=1,N=1lat X =0, 
g¢=1,0=0,N=1lat X =1, 
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Local solid Nusselt number, Sherwood number of nano particles and Sher- 
wood microorganism are defined as 


Nuy = —(28)x=0,1, Shy = —(28)x~0,1, and Nny = —($%)x-=0,1- 
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Figure 1: Physical geometry and coordinate system. 


The average quantities of Nusselt number, nanoparticle Sherwood num- 
ber, and microorganism Sherwood number is defined as 

Nuavg =fy Nuy dY, 

Shavg =Jy Shy dY, 

Nnavg =Jy Nny dY. 


3 Numerical method 


To find the numerical solution to (8)—(11), partial differential equations with 
boundary conditions Galerkin’s weighted residuals finite element method are 
solved with the help of MATLAB [9]. In this method, a two-dimensional 
field is divided into small triangular parts, in which each part is named an 
element. Over each element, assume a piecewise trial function. 

Let U, 6, ¢, and N be approximated by UV = S7?_, UiE;, 0 = 732, Hi, 
o= 5 o&, and N = 4 Ni€i, where €; is the linear interpolating 
functions over each triangular element. The FEM model matrix is as follows: 
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- ee the system of equations, the functions are incorporated, which 
are assumed to be known. After applying the boundary conditions, a matrix 
of system of linear equations is formed, which is solved by using the Gauss— 
Seidel iteration method. The convergence of the solution is assumed when the 
relative error for each variable between two consecutive iterations is observed 
below the convergence criteria such that |y"*t — "| < 107°, where n is the 
number of iterations and w stands for V,0,C. 


To choose the grid size, the grid independence test is performed for 
21 x 21,41 x 41,61 x 61,71 x 71,81 x 81,91 x 91 grid sizes. The grid in- 
dependence test reveals that the grid size 81 x 81 is sufficient to study in the 
the bioconvection phenomena. 


Table 1: Accuracy test of mean Nusselt numbers with the literature. 
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Authors Ra=10 | Ra = 100 | Ra = 1000 
Varol, Oztop, and Pop [33] 13.564 
Cross, Bear, and Hickox [12] - - 13.470 
Manole [22] - 3.118 13.637 
Sheremet and Pop [28] 1.079 3.115 13.667 
Present results 1.081 3.1271 13.715 
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4 Result and discussion 


The present equations (8)—(11) are numerically investigated and analyzed in 
the porous square cavity filled with nanofluid and gyrotactic microorganisms 
at different inclination angles. The numerical investigation is carried out with 
the following parameters considered Rayleigh number (Ra = 25), bioconvec- 
tion Rayleigh number (Rb = 15), Lewis number (Le = 1), Peclet number, 
thermophoresis parameter, Brownian parameter, buoyancy ratio parameter, 
Schmidt number 0.1, Prandtl number 6.9, and inclination angle (6 = 30°) 
unless when it is mentioned. Streamlines are presented in Figure 2 for var- 
ious values of inclination angle (0). When 6 = 0°, the cells moved in the 
clockwise direction in the square enclosure. The absolute maximum of the 
stream function is |Wyax| = 1.56. At 30°, inclination angle the maximum 
stream function value is |Wj,ax| = 2.0435, and the rotation of the cell moves 
in the same direction. The inclination angle increased to 90°, and the veloc- 
ity of the fluid flow is reduced. At 120°, the intensity of the flow increased to 
|Wnax| = 2.0421. In this, the cell moves to the center regime with the inten- 
sity of the gravitational force. At 150° inclination angle, the flow strength 
is low and the cell moves near to corners of the left down wall and right 
top walls. At 180°, the fluid flow velocity is less with the strength of cell 
|Wnax| = 1.5629. 

Isotherms are demonstrated in Figure 3 for various inclination angles from 
0° to 180°. At inclination 30°, the temperature distribution is indicating the 
stratified diagonally. At 180°, isotherms are parallel to the vertical walls, 
which shows the transfer of heat in the mode of heat conduction. 

Nanoparticle isoconcentrations of the fluid for various angles are expressed 
in Figure 4. At 0° to 30° inclination angle, the nanoparticle volume fraction 
was raised near the bottom cavity wall. When the square enclosure was 
moved from angle 90° to 120°, the nanoparticle isoconcentration decreased. 
When the angle is 120° to 150°, again the concentration of nanoparticles 
increases. When the angle of inclination is inclined from 150° to 180°, the 
concentration of nanoparticles is decreased. In these all angles, the cell is 
divided into two different parts: one is formed near the bottom adiabatic 
wall and the other part is formed as a semi-opened vertex close to the top 
adiabatic wall. 

Figure 5 displays isoconcentrations of microorgnisms for the various incli- 
nation angles. At 0° to 30° inclination angle, two types of cells are formed: 
one cell is near the bottom adiabatic wall and the cell moves from the cold 
wall to the hot wall, the second is an open semi vertex formed at the top 
cavity wall, and it moves from the hot wall to the cold wall. When the 
inclination varies from 0° to 180°, the movement of cells in the isoconcentra- 
tions of microorganisms and nanoparticle volume fraction is the same. The 
concentration of microorganisms’ maximum value is found at 6 = 30°, 120°. 

In Figure 6, the effect of Rb and Nt on average Nusselt number, Sher- 
wood numbers of nanoparticles, and Sherwood number of microorganisms is 
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Figure 2: Streamlines for the inclination angle 6 = 0° — 180° 


discussed. Moreover, Rb increases Nuaug and Shavg from 0° < 6 < 90°, 90° 
< 6 < 180°, and it reaches high. Indeed, at 0°, 90°, and 180° the average 
Nusselt number and average Sherwood number values are low. Also, Rb in- 
creases Nnavg from 0° to 180°. In addition, Nt increases Nuayg from 0° to 
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Figure 3: Isotherms for the inclination angle 6 = 0° — 180° 


180°, and it shows wavy behavior, and Shavg and Nnavg are also increased 
with the increase of thermophoresis parameter. 
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Figure 4: Isoconcentration of nanoparticle volume fraction for the inclination angle 
6 = 0° — 180° 


5 Conclusion 


The effect of porous square cavity inclination with the horizontal surface with 
nanofluid and gyrotactic microorganisms was analyzed with the streamlines, 
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Figure 5: Microorganism isoconcentrations for the inclination angle 6 = 0° — 180° 


isotherms, nanoparticle volume fraction, and microorganism isoconcentration 
from 0° to180°. 


1. The velocity of the nanofluid flow is high at the angle 30°, 120° and in 
the remaining angles, the flow intensity is low. 
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Figure 6: Representation of (a) Average Nusselt number (b) Average nanoparticle Sher- 
wood number (c) Average Microorganism Sherwood number for angle versus Rb and 
angle versus Nt 


2. The temperature distribution of the nanofluid is affected by the square 
cavity inclination. 

3. Nanoparticle isoconcentration and microorganism isoconcentrations 
are high at the 30° and 120°, and at the remaining angles, the value is low. 
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4. Thermophoresis parameter increases Nttayg, Shavg, and Nrayg from 


0° <6 < 180°. Also, Nt increases Nitgyg from 0° < 6 < 180°, but at 0° and 
180°, the value is low. 5. Bioconvection Rayleigh number increases Nuayg, 
Shavg, and Nnavg- 
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